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ABSTRACT 


A growing avenue for determining the prevalence of life beyond Earth is to search for “technosig- 
natures” from extraterrestrial intelligences/agents. Technosignatures require significant energy to be 
visible across interstellar space and thus intentional signals might be concentrated in frequency, in time, 
or in space, to be found in mutually obvious places. Therefore, it could be advantageous to search 
for technosignatures in parts of parameter space that are mutually-derivable to an observer on Earth 
and a distant transmitter. In this work, we used the L-band (1.1-1.9 GHz) receiver on the Robert C. 
Byrd Green Bank Telescope (GBT) to perform the first technosignature search pre-synchronized with 
exoplanet transits, covering 12 Kepler systems. We used the Breakthrough Listen turboSETI pipeline 
to flag narrowband hits (~3 Hz) using a maximum drift rate of +614.4 Hz/s and a signal-to-noise 
threshold of 5 — the pipeline returned ~ 3.4 x 10° apparently-localized features. Visual inspection 
by a team of citizen scientists ruled out 99.6% of them. Further analysis found 2 signals-of-interest 


Corresponding author: Sofia Sheikh 
ssheikh Qseti.org 


that warrant follow-up, but no technosignatures. If the signals-of-interest are not re-detected in future 
work, it will imply that the 12 targets in the search are not producing transit-aligned signals from 1.1 
— 1.9 GHz with transmitter powers >60 times that of the former Arecibo radar. This search debuts a 
range of innovative technosignature techniques: citizen science vetting of potential signals-of-interest, 
a sensitivity-aware search out to extremely high drift rates, a more flexible method of analyzing on-off 
cadences, and an extremely low signal-to-noise threshold. 


Keywords: technosignatures — search for extraterrestrial intelligence — biosignatures — astrobiology 


1. INTRODUCTION 


A vast and still-growing part of our astronomical exploration is the search for life elsewhere in the universe. Many 
programs look for such life on exoplanets through their biosignatures, surface features (e.g., Coelho et al. 2022), 
or atmospheric constituents (e.g., Thompson et al. 2022) that indicate the presence of biological activity. Another 
complementary strategy is to look for the technosignatures of technologically-capable life — Extraterrestrial Agents 
(ETA)s, to use the terminology of Dobler & Raab (2021) — which may be more abundant, long-lived, highly-detectable, 
and unambiguous than other previously-described biosignatures (Wright et al. 2022). 

The most popular technosignature search strategy to date is radio searches for artificial emission (as pioneered by 
Drake 1961), which has grown exponentially in recent years, making use of cutting-edge computational techniques 
and newly-developed hardware (e.g., Harp et al. 2018; Ma et al. 2022). However, even with the renewed observational 
energy, the search space remains mostly unexplored (Wright et al. 2018a). This provides an opportunity for radio 
observation projects, large and small, to make an impact by filling in unexplored regions of parameter space. 

One suggestion on how to best navigate this huge parameter space is to use “Schelling Points” (Schelling 1960; 
Wright 2020), to prioritize mutually-derivable parts of parameter space which a transmitter and receiver can both 
land upon without any prior communication. This allows for more efficient traversal of parameter space — potentially 
leading to a technosignature detection much sooner — and also can be more power efficient for the transmitter, which 
can focus its energy in particular directions and times. One application of this idea is to prioritize certain places and 
times synchronized with an astronomical event (Pace & Walker 1975). In some of these synchronized strategies, the 
event is external to the transmitter and receiver’s systems, i.e., a nearby supernova or nova (Tang 1976; Lemarchand 
1994; Davenport et al. 2022). In another application, the synchronizing event is some feature of the transmitter or 
receiver’s system, observable or predictable by both parties (e.g., Corbet 2003). 

Here, we make use of exoplanetary transits as temporal Schelling points (Kipping & Teachey 2016): if an ETA 
transmitter on an exoplanet sends a signal towards its anti-stellar point, the signal will necessarily arrive to any observer 
along the line-of-sight at the same time as the exoplanet appears to be at mid-transit. This provides a known time and 
place (during an exoplanetary transit) for the observer to look for a signal. The transmitter may be targeting a specific 
system that happens to lie within its ecliptic, and thus sends a signal once each of its exoplanetary years. Conversely, 
the transmitter may be constantly targeting its exoplanet’s anti-stellar point, sending out a transmission which sweeps 
across its ecliptic. In the extreme case, a tidally-locked planet hosting a single surface-locked transmitter could be 
constantly sweeping the ecliptic via the anti-stellar point, perhaps powered by energy collected photo-voltaically from 
the star-facing side of the exoplanet. Observations at this special periodic epoch offer a way to sample the large 
possibilities of repetition rates associated with periodic transmissions (Wright et al. 2018b). Completely divorced from 
this Schelling Point logic, ETAs may preferentially emit high-power microwave radiation in their ecliptics (as we do in 
our solar system due to space communications), making transiting systems potentially more favourable for detecting 
unintentional artificial signals. 

In this work, we perform the first radio technosignature search of exoplanet-hosting stars, where observations were 
pre-synchronized with their transits. This is a complementary approach to that performed by Franz et al. (2022), 
which looked back into archival data to identify serendipitous observations during transit. We also follow the growing 
tradition of conducting technosignature work and obtaining novel astronomical results using a cohort-based, academic 
course centered research model (Margot et al. 2021; Tusay et al. 2022). 

In Section 2, we describe how the targets were chosen. In Section 3 we discuss the observing plan, the observation 
parameters, and the data format. We cover the methods for the search, including the software, assumptions, and 
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Figure 1. Confirmed transiting exoplanets in the Kepler field as observable from the Green Bank Telescope (GBT) during on 
March 25, 2018. The horizontal errorbars depict the transit duration, whereas the colour of the markers represent the stellar 
effective temperature. Points without visible errorbars have extremely well-defined transit midpoints. The grey region marks 
our observing window, whereas the targets observed are circled in black. We did not include a habitability requirement in planet 
selection. 


chosen search parameters in Section 4. We discuss our findings and upper limits in Section 5. Finally, we discuss and 
conclude in Section 6 and Section 7 respectively. 


2. TARGET SELECTION 


We first compiled a list of 540 exoplanets that transit during a 2-day potential observing window (March 23, 2018 
— March 25, 2018) as a function of distance, using ephemerides from the NASA Exoplanet Archive (Figure 1; Akeson 
et al. 2013). Down-sampling from this subset of transiting exoplanets, we selected 436 planets discovered by the Kepler 
mission due to their limited range of celestial coordinates. This enabled us to minimize our slew times during the 
nodding sequence, and boost the observing duty cycle. While there have been previous searches focused on transiting 
exoplanets from the Kepler mission (Harp et al. 2016; Siemion et al. 2013b), these have not prioritized observations 
during planetary transits. As explained above, we consider the transit a temporal Schelling point and hence aim to 
maximize the fraction of transit that we observed. To further improve our observing efficiency, we decided to observe 
alternate transiting planets as part of the on-off-on sequence used to identify Radio Frequency Interference (RFI), and 
identified pairs of exoplanets in the Kepler field that are transiting at similar times and are particularly close on-sky!'. 
Of the Kepler transiting planets in our observing window we down-selected to 12 which are a) fairly close in their 
positions on-sky to allow for efficient alternating (on-off) exposures, and b) would result in sufficient sampling both 
in-and-out of transit, to fill the 6 hours of accepted observing time from Breakthrough Listen (BL).. The properties 
of our 12 observed targets are listed in Table 1, and their on-sky positions are visualized in Figure 2. 


3. OBSERVATIONS AND DATA 


We pre-planned our observing sequence to ensure that we would hit each exoplanet during its transit, using individual 
5-minute observing blocks spanning our 6 hour observing window on March 25, 2018 (starting at 11:00 UT); we elected 


l None of these pairs are close enough on the sky to potentially cause source confusion, i.e., the angular distance between the targets >> 
GBT beam diameter in L band. 
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Figure 2. A plot of all exoplanetary systems (blue stars) observed with the GBT for this project, overlaid on the Kepler field 
of view, which covers 115 square degrees on sky (Mullally et al. 2016). The name of each system as it appears in Table 4 is 
shown. Boyajian’s star (an off-source target) is also included for reference. 


to use the standard BL Green Bank Telescope (GBT) integration of 5 minutes, and assumed a 2 minute overhead 
to account for e.g., the slewing and settling of the telescope. Table 4 shows the log of observations, including their 
relative temporal position to the target exoplanet’s mid-transit point. We created individual GBT observing scripts 
for each pair which toggled back-and-forth between the two targets until they were replaced by the next target script. 


Table 1. Stellar and planetary properties for the twelve transiting Kepler planets observed in this work. We use the properties 
from Morton et al. (2016) for all planets, except Kepler-446b which we pull from Muirhead et al. (2015). 


Target RA Dec Distance Radius Period to Teg Teff 

(J2000) (J2000) (pc) (Re) (da) (BJD - 2454000) (K) (K) 
Kepler-446b 18h49m0.05s 44d55m15.96s 120 1.5 +0.25 1.57 965.91 648 3464 
Kepler-537b 19h19m31.20s 51d16m48.00s 465 1.41499% 3.25 1004.28 1181 5703 
Kepler-723b 18h59m19.32s 44d39m29.30s 965 12.194538 4.08 1002.64 1016 5655 
Kepler-732c 18h54m55.68s 45d57m31.57s 150 Lore 0.89 967.54 893 3582 
Kepler-738b 19h10m16.73s 46d34m4.30s 992 apres. 24.59 1006.70 530 5474 
Kepler-842b 19h29m16.80s 49d38m60.00s 552 16". 1.22 966.46 1252 4838 
Kepler-992b 18h57m43.20s 47d38m24.00s 268 1.62+9:94 20.16 977.34 510 4975 
Kepler-1039b 19h53m16.80s 45d18m36.00s 324 1.464993 0.93 964.58 2080 4777 
Kepler-1053b 19h25m40.80s 39d7m48.00s 171 0.981903 2.41 965.43 953 4507 
Kepler-1164b 19h10m7.20s 38d53m24.00s 447 1,121007 -3.98 966.62 936 5143 
Kepler-1222b 19h39m33.60s 49d22m48.00s 455 0.79 0.06 1.92 965.47 1246 5083 
Kepler-1332b 19h24m7.20s 43d54m36.00s 465 1.371999- 11:87 973.21 728 5523 
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This sequence was adjusted dynamically throughout the observing session to select targets closest to mid-transit, in 
light of actual slew times and unanticipated observing overheads. For example, we started our transit observations 
with Kepler-992b for scans 0010 and 0012, the former of which spanned the transit midpoint. During these scans, we 
observed Kepler-738b as our ‘off’ target. After covering the transit midpoint for Kepler-992b, we switched to Kepler- 
1039b and Kepler-732c, which had the next occurring transit midpoints. The total in-transit time for each planet in 
the sample varied in duration from 0.65-3.9 hours (median: 1.8 hours), so the 5-minute scans covered approximately 
5% of each transit. 

We followed a similar logic through the rest of our observing window, where we tried to observe targets for at 
least 2 scans each during transit; furthermore, we prioritized observations of transit midpoints and tried to minimize 
slew times. We bracketed our observing block with calibration observations of quasar 3C 295 (scans 0006 and 0007) 
in the beginning, and pulsar B0329+54 (scan 0059) at the end. Additionally, we obtained one scan (scan 0009) of 
KIC 8462852, commonly known as Boyajian’s star (Boyajian et al. 2018) before starting our transit sequence; it was a 
conveniently-located off-source also targeted by BL laser technosignature searches using Lick Observatory’s Automated 
Planet Finder telescope (Lipman et al. 2019). 

Data were recorded using the L-band receiver (1.1-1.9 GHz) (Balser et al. 2021) and the BL backend on the GBT, 
which in 2018? was capable of digitizing up to 6 GHz of instantaneous bandwidth in 2 polarizations at 8-bit resolution 
— for more information, see MacMahon et al. (2018). The raw voltage data is then channelized using the gpuspec 
code to produce three spectral SIGPROC filterbank files containing Stokes I total power /frequency/time datacubes, each 
file at a different time and spectral resolution (Lebofsky et al. 2019). In the following analyses, we make use of the 
high spectral resolution data product, with a frequency resolution of ~ 3 Hz and a time resolution of ~ 18s. 

We performed data quality checks using the calibrator observations at the beginning and end of the observing session. 
The pulsar B0329+54 was easily visible in a prepfold plot (Ransom 2011), providing a first-order confirmation that 
the system was working as expected. In addition, it had an expected Signal-to-Noise Ratio (S/N) of 9172 given our 
system parameters, observing parameters, and its characteristics in the ATNF pulsar database, and was detected at 
an S/N of 5306 (57.85%). Given that pulsars in general (and this pulsar in particular) are not perfect flux calibrators 
(they exhibit variability in flux over time due to scintillation, which was distinctly present in the scan), this is entirely 
within the range of expected outcomes. 

We also used the first 3C 295 scan as a continuum flux calibrator “on”-source, and the following observation of 
Boyajian’s Star as an “off” source, to derive a system temperature. Assuming a flux density of 22.22 Jy at 1408 MHz 
(as measured by Ott et al. 1994), and a spectral index of a = —0.7, as found for the calibrator’s hotspots in a 
recent LOFAR observation (Bonnassieux et al. 2022), we obtain an empirical system temperature measurement of 
Tsys = 22.66 K. The theoretical value given in the GBT Proposer’s Guide for the L-band receiver is Tsys ~ 20 K, which 
is consistent with the experimental results. 

The data from this campaign is hosted at UC Berkeley, managed by the BL project, and globally accessible via http 
at https://bldata.berkeley.edu/kepler-transits/. The repository includes the high-resolution spectral data product 
used for the following analysis, as well as the spectral line-resolution, and high-resolution data products from the same 
observations. Specifications for these other data products can be found in Lebofsky et al. (2019). 


4. SEARCH METHODS AND DATA REDUCTION 


BL has a well-established data pipeline for performing narrowband Search for Extraterrestrial Intelligence (SETI) 
searches on high frequency resolution filterbank files. This pipeline consists of import, plotting, and other utility 
functions in the blimpy package (Breakthrough Listen Collaboration 2019; Price et al. 2019), and a narrowband 
search code turboSETI (Enriquez & Price 2019) based on the de-dispersion algorithm of Taylor (1974) and adapted 
for de-Doppler searches. These pipelines have frequently been used for SETI searches in the recent literature (e.g., 
Siemion et al. 2013a; Smith et al. 2021). 

For this work, we ran turboSETI’s narrowband hit-finder with a maximum drift rate of + 614.4Hz/s and a S/N 
threshold of 5. Both of these parameters are unusual choices, for the reasons presented below. 

A typical drift rate for radio SETI searches is ~ +10 Hz/s, within a factor of a few (e.g., Sheikh et al. 2020; Franz 
et al. 2022). However, Sheikh et al. (2019) showed that known exoplanetary systems could produce drift rates up to 
+ 200 nHz (i.e., +200 x v Hz/s, where v is the observing frequency in GHz). The largest drift rates would be expected 


2 It is now capable of digitizing up to 12 GHz of instantaneous bandwidth. 
3 https: //www.atnf.csiro.au/research/pulsar/psrcat/ 
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Figure 3. A histogram showing the frequency distribution of the hits (lighter) and events (darker) in a logarithmic (left) and 
a linear (right) scale. The notch filter is responsible for the absence of hits between 1200-1340 MHz. Three distinct frequency 
regions contain more than 100000 hits each: 1165-1185 MHz, 1375-1385 MHz, and 1570-1580 MHz. These three regions are 
discussed in further detail in Section 5.1. 


from exoplanets with the largest radial accelerations relative to the receiver on Earth, e.g., transiting exoplanets that 
are viewed edge-on. This work marks the first time that a SETI survey has followed the 200nHz recommendation, 
by using a maximum drift rate of +614.4Hz/s; this drift rate is sufficient to capture drift rates of +200 nHz even at 
1.9 GHz, the highest frequency in these observations. 

It should be noted that turboSETI, in its current configuration, does not automatically account for the sensitivity 
loss sustained when incoherently-dechirping past the “one-to-one point” Vi—to—1, where a signal drifts one frequency 
bin in every time bin (Sheikh et al. 2019). Here, we implement the first of the two partial solutions to this problem 
described by Margot et al. (2021). We search the original high spectral-resolution filterbank file with a drift rate range 
from 0 to +11_to-1 Hz/s, in steps of ~ 0.01 Hz/s. Then, we “scrunch” the file in frequency by a factor of 2, halving the 
number of bins by summing every other bin with its following neighbor. We then search again, using a new drift rate 
range of +1;_to0-1 to +2 x Vi—to—1 Hz/s (the new one-to-one point in the new file), with a correspondingly doubled 
drift rate step. We repeat this process until we have covered the entire desired drift rate range, which requires a series 
of 12 scrunch-and-search iterations. It should be noted that each frequency-scrunch, though it maximizes the S/N 
within its range of drift rates, still causes an unavoidable V2 loss in sensitivity due to the original dimensions of the 
incoherent sum that produced the data product. The sensitivity losses are discussed further in Section 5.4, and it 
should be noted that exact signal positions and drift rates within the band may also cause irregularity in sensitivity 
with this method. 

Recent SETI searches have often chosen an S/N of 10 for their hit thresholds, including Price et al. (2020) and 
(Margot et al. 2021). In narrowband radio signal detection these limits are primarily dictated by the filtering resources 
of the search, rather than any inherent statistical significance — the environment is so contaminated by RFI that false 
alarm rates from other common astronomical distributions, e.g., white noise or Gaussian statistics, do not apply. In 
this work, we instead decided upon a lower S/N of 5, which allows us to double the sensitivity of the search. It should 
be noted that this causes an immense number of hits to pass the turboSETI filtering step, in addition to the shorter 
BAB cadence described in the following section: we managed this step with citizen science (Section 5.2), but also note 
its difficulties (Section 6). 


5. RESULTS 


We ran turboSETI’s hit-finding algorithm on every observation in Table 4. We ignored the region from 1.20- 
1.34 GHz, which corresponds to the L-band “notch filter”, which applies to the most RFI-contaminated region of the 
spectrum at the GBT. This generated a total of ~2.53 million hits. We then used turboSETI’s event-finding capability 
to compare hits in each Kepler observation (and Boyajian’s star) to hits in each of the observations directly preceding 
and following it; this is the equivalent of an off-on-off or BAB cadence. In rare cases where the same target was 
observed twice in a row to realign the observations with the transit timeline (e.g., scans 0018 and 0019, both targeting 
Kepler-738b), we used the next closest scan in the preceding/following direction that was on a different target. This 
event-finding process resulted in 338473 unique events. The frequency distribution of the hits and events are shown 
in Figure 3. 
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Figure 4. a) An example of an event in the 1165-1185 MHz region, consistent with GPS L5. This event is RFI, as it appears 
in multiple different targets in consecutive observations. The dashed red line represents turboSETI’s best-fit drift rate for the 
detected hit. b) An example of an event in the 1375-1385 MHz region. These events occur for only a fraction of the observation 
(here, about 100 seconds) within a single scan — this makes it difficult to determine whether they are impulsive RFI or whether 
they are true transient signals localized on the sky. However, because we observe the same morphology of signal in multiple 
targets, and because of the match to the GPS L3 downlink, we can assign these events as RFI. c) An example of an event in 
the 1570-1580 MHz region. These events are similar to the example shown in subfigure b in that they are degenerate between 
transient and localized signals; this is a common challenge for single-dish technosignature searches. The same signals were 
identified in multiple targets, however, indicating that they are indeed RFI — likely the GPS L1 downlink. d) An example of a 
high-drift rate signal around 1626 MHz. We attribute this event to an Iridium satellite. It is clearly present in all three panels, 
and therefore RFI, but the changing slope of the signals over time confounds the linear search algorithm. 


5.1. RFI-Heavy Bands 


We find that three frequency ranges — 1165-1185 MHz, 1375-1385 MHz, and 1570-1580 MHz — contain the majority 
of the RFI in our observations: 5% of the band (excluding notch filter) accounts for 56% of the hits. These ranges are 
consistent with the interference-heavy regions discussed by Price et al. (2020). Here, we briefly discuss each of these 
frequency ranges in turn, and attempt to associate them with Federal Communications Commission (FCC) frequency 
allocations’. 

1165-1185 MHz: There is an FCC frequency allocation for aeronautical radionavigation and radionavigation- 
satellite (including space-to-Earth) transmissions between 1164-1215 MHz, covering this observed interference region 
in Figure 3. The distinct peak from 1165-1185 MHz is consistent with the GPS L5 downlink®. An example of an event 


in this region that passed turboSETI’s thresholds is shown in Figure 4a. 


4 https: //transition.fcc.gov /oet /spectrum/table/fectable.pdf 


5 https: //www.nist.gov /pml/time-and-frequency- division /popular-links/time-frequency-z/time-and-frequency-z-g¢ 


1375-1385 MHz: This narrow region of the spectrum is occupied by the GPS L3 downlink band, which provides 
communications links for the U.S. Nuclear Detonation Detection System. An example of an event in this region that 
passed turboSETI’s thresholds is shown in Figure 4b. 

1570-1580 MHz: Once again, this interferer falls within an FCC allocation dedicated to aeronautical radionavi- 
gation and radionavigation-satellite (including space-to-Earth) transmissions; in this case, however, the allocation is 
much wider (1559 MHz-1610 MHz) than the region where we see the majority of interference. An example of an event 
in this region that passed turboSETI’s thresholds is shown in Figure 4c. Given the frequency range, and the presence 
of other GPS downlinks within the dataset, these hits can likely be attributed to the GPS L1 downlink centered at 
1575 MHz. 

Finally, we detected a series of swooping signals (changing from high drift rate, to low drift rate, to high again) 
between 1616—1626.5 MHz that also passed turboSETI’s event filtering. We highlight these signals due to their high 
drift rates, which reached an absolute value of nearly 70 Hz/s.These signals account for the event spike above 1600 MHz 
in the right panel of Figure 3. An example is shown in Figure 4d. We attribute these signals to the Iridium satellite 
constellation’s L-band downlink®, and note that these signals are bright enough to have passed turboSETI event filters 
in multiple GBT L-band SETI observations, even those with lower drift rate bounds. (e.g., Enriquez et al. 2017; 
Tusay et al. 2022). 


5.2. Filtering By Citizen Scientists 


At this stage, with > 340,000 events, most radio technosignature campaigns would make a change to the event 
thresholds to get them to a number suitable for visual inspection by a single researcher. For example, the S/N could 
be raised from 5 to 10 (leaving only 92,000 events), or even 25 (leaving only 26, 000 events), or the drift rate could be 
further reduced. However, making these changes would lower our sensitivity. For this campaign, we decided instead 
to apply the power of crowdsourcing via a limited citizen science project. 

We created .pdf files containing 1000 plots each, making them accessible via cloud services such as Box and Google 
Drive, and developed an hour-long lecture and interactive quiz to train volunteers. Six undergraduate volunteers from 
the Penn State Pulsar Search Collaboratory (PSPSC) club” looked through output plots such as those shown in Figure 
4, and flagged any signal that appeared in only the on-source middle panel. Due to the COVID-19 pandemic, the 
PSPSC members only completed approximately 20% of the sample in the Fall 2020 semester — additional volunteers 
were recruited from the lead author’s professional network to assist with signal filtering. 

In the end, 13 citizen scientists (named in the Acknowledgements) flagged approximately 0.4% (~ 1600 signals) of 
the dataset for further analysis, reducing the number of interesting events to a few thousand. It should be noted that 
this number is approximate: some signals were identified multiple times at different drift rates (see Section 6), and 
clusters of plots with repeating signals (such as those in Figure 4a) were grouped as a single phenomenon. Due to 
the size of the dataset compared to the number of volunteers, we did not have the resources to send the same data 
to multiple recipients, as is done in projects such as Zooniverse®. To ensure quality in flagging, the majority of the 
work was done in collaborative “hack sessions” with multiple participants on the same call, so plots-of-interest or 
borderline cases could be discussed and viewed by multiple participants simultaneously. The citizen scientists’ effort 
meant that the filtered dataset was now at a manageable size for visual inspection by the authors, all while maintaining 
the extremely sensitive S/N 5 threshold. 


5.3. Filtering Events with Context 


Each of these plots, in isolation, had the characteristics of a localized signal on the sky (only present in the on-source 
observation). However, following the technosignature verification framework of Sheikh et al. (2021), these signals 
should also be analyzed in the context of similar, nearby signals. We grouped the remaining few thousand plots by 
frequency and morphology. Then, we checked to see if identical signals appear in two different targets (indicating that 
they are not local to either) or if we see a near-identical signal within a set that is definitely interference (the way that 
blc1 was eventually disproved by Sheikh et al. 2021). With this strategy, we reduced the pool of signals-of-interest to 
479. The frequencies, drift rates, and S/Ns of these 479 signals are shown in the context of the overall hits and events 
in Figures 5, 6, and 7 respectively. 


6 https: //apollosat.com /iridium-satellite-frequency- bands 
T An affiliate of the Pulsar Search Collaboratory project described by Blumer et al. (2020) and others. 


8 https: //www.zooniverse.org/ 
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Figure 5. A histogram showing the frequency distribution of the hits (light gray), events (light purple), and signals-of-interest 
(purple). The signals-of-interest are distributed relatively evenly throughout the spectrum. Note the gap from 1.20-1.34 GHz 
due to the GBT L-band receiver’s notch filter. 
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Figure 6. A histogram showing the drift rate distribution of the hits (light gray), events (light orange), and signals-of-interest 
(orange). The signals-of-interest are found only at low absolute drift rates. 


10 


"~~ Events 
10° EEE Signals-of-Interest 


104 


Count 


107 
10? 
10? 


10° 


0 2000 4000 6000 8000 10000 
Signal-to-Noise Ratio 


Figure 7. A histogram showing the S/N distribution of the hits (light gray), events (light green), and signals-of-interest 
(green). The signals-of-interest are primarily faint. There are also 5 hits with S/N > 10000 which are not shown on this plot 
for readability. 


It is, however, possible that the 479 signals-of-interest have disqualifying features elsewhere in the observing session; 
for example, a sub-threshold repetition of the signal on a non-adjacent scan would not be picked up in the previous 
filtering step. Therefore, as a final verification, we plotted each of the remaining 479 plots in the context of the entire 
session of observations (all 52 observations in Table 4), to produce stacked waterfall plots with extreme aspect ratios. 
We eliminate a signal-of-interest from our list if it is a continuation or repetition of signals in other targets at any time 
during the morning of observations: this is a similar strategy to that applied in Section 5.2 but now across the largest 
possible time baseline. 

This process eliminated all but 2 signals-of-interest: one in a scan of Kepler-1332b and one in a scan of Kepler-842b. 
Neither of these signals are actually narrowband, but they were bright enough to be detected by turboSETI and we 
did not want to restrict ourselves to only the signal morphologies we expected, if we serendipitously found a signal 
with a different morphology. Interestingly, both the detection in Kepler-1332b (during scan 0031) and the detection 
in Kepler-842b (during scan 0056) were during their respective transit midpoints. The waterfall plots are shown in 
Figures 8 and 9, and the signal properties are summarized in Table 2. 


Target Scan | MJD of Scan Frequency Transit Phase Signal Type 
Kepler-1332b | 0031 | 58202.59350 | 1749.4209 MHz | Transit Midpoint Pulse 
Kepler-842b | 0056 | 58202.69226 | 1040-1438 MHz | Transit Midpoint | Broadband Pulse 


Table 2. Properties of the two interesting signals-of-interest reported in this study. “Transit Midpoint” indicates that the 
midpoint of the planet’s transit occurred during this scan. 


As both signals-of-interest are contained within a single time bin in the high frequency-resolution product, we then 
checked the mid-resolution and high time-resolution products for evidence of the same signal. The high time-resolution 
did not reveal anything of interest, likely due to its low frequency-resolution, but the mid-resolution product (~3 kHz 
frequency resolution, 1 second time resolution) does provide additional information. The signal-of-interest in Kepler- 
1332b splits into two still-unresolved bursts spaced about one second apart. The signal-of-interest in Kepler-842b is 
also revealed to have a more complex structure in time, also periodic in frequency, which repeats across the entire 
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Figure 8. The signal-of-interest detected in Kepler-1332b. 


bandwidth. These mid-resolution plots are shown in the Appendix. This new context about the morphology of the 
signal-of-interest lends more evidence to the idea that the signal-of-interest is broadband RFI. 

We are close to the limit of what we can deduce from this dataset for these two signals-of-interest. With the 
information we have, it is impossible to exclude these signals as being due to RFI, which should be treated as the most 
likely explanation.. Therefore, we recommend a follow-up observation of these two systems in L-band, during transit, 
with a different instrument. We are planning to do this with the newly-refurbished Allen Telescope Array (ATA) (e.g., 
Farah et al. 2021). These signals-of-interest do not merit the scrutiny of blc1 (Sheikh et al. 2021) because they only 
appear in a single 18 second integration; while this may be expected for synchronized transmitters which only send a 
signal at mid-transit, it also makes the standards-of-evidence higher (no proof of signal localized on sky, no drift rate 
measurement, etc.), requiring reobservation as the next step before more detailed analyses. 


5.4. Upper Limits 


We calculate technosignature upper limits following the method of Price et al. (2020). The minimum detectable flux 
density, Fmin for a narrowband signal is given by: 


2kg Teys ôv 
a Aeff Npoltobs 


(1) 


Here, S/N nin is the signal-to-noise threshold, kg is the Boltzmann constant, Tsys is the system temperature, Aeff is 
the effective collecting area of the telescope, dv is the frequency channel width, and npo is the number of polarizations. 
We make the same assumption as Price et al. (2020): that the transmitting signal is the same width as the frequency 
channel dv, so no additional constant factor is needed to downweight Fmin- 

With an S/N min of 5, a Tsys of 20K, an Aer of 5655 m? (equivalent to using an effective diameter of 100m for the 
GBT, with an L-band aperture efficiency of 0.72), a dv of 3 Hz, an npol of 2 and a tops of 300 seconds, we calculate 
an Finin of 3.45 x 107~?° W/m?, equivalent to 1.15 Jy. This is consistent with other BL searches at L-band with the 
GBT (e.g., Enriquez et al. 2017; Price et al. 2020), but with extra sensitivity due to the low signal-to-noise threshold. 
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Figure 9. A portion of the signal-of-interest detected in Kepler-842b. This signal-of-interest spans 1040-1438 MHz, so this is 
not the full signal, but rather a narrow bandwidth example of one of the “hits” on this signal generated by turboSETI. 


The population of Kepler exoplanets is relatively distant, with all of the targets in this study falling between 150- 
992 pc. The Equivalent Isotropic Radiated Power (EIRP) is calculated as: 


EIRPpin = 4rd? Frin (2) 


where d is the distance to the furthest target in the survey. Using 992 pc, we get a survey-wide EIRP min of 406 TW, 
or ~ 20x the EIRP of the Arecibo radar. The two targets that showed signals-of-interest, Kepler-1332b and Kepler- 
842b, are located closer, at 465 pc and 552 pc respectively. This leads to target-specific EIRPs of 89 TW and 125 TW, 
or 4.5-6.3x the EIRP of Arecibo. 

Finally, one can also express an upper limit in the form of a fraction of the 8-dimensional “haystack”, as described 
by Wright et al. (2018a). Using the same minimum and maximum axis bounds as described in that work, we calculate 
a total haystack fraction of 4.42 x 107° covered by this search, comparable to Project Phoenix searches with Parkes 
and Arecibo (e.g., Backus et al. 2002). 


6. DISCUSSION 


This search represents the lowest S/N ever chosen for a Breakthrough Listen affiliated project on the GBT. However, 
caution is advised if following our procedure. We found that lowering the S/N to 5 made turboSETI’s temporal on-off 
filtering much less effective. The lower S/N magnified the likelihood of e.g., a bright pixel in the “on” target causing a 
hit to be flagged as an event. It may be more effective to recoup sensitivity in other parts of the signal processing chain, 
rather than as a threshold in turboSETI. This may be especially relevant for high drift rate signals (using methods 
such as boxcar convolution, e.g., Adámek & Armour 2020). 

In addition, vetting ~ 3.4.x 10° events involved a large amount of work, even considering the citizen science approach. 
Using a single on-source with two off-sources produces a degeneracy between intermittent RFI and sky-localized signals 
in a single-dish telescope, compounding the issue. We recommend that future searches invest in a better understanding 
of their RFI environment and try to algorithmically filter further before performing visual classification. 
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One approach is to better characterize the frequencies, drift rates, and morphologies of common interference sources, 
and only down-select signals that are not particularly consistent with those properties. Another approach is to use 
employ spatial filtering with a multi-beam receiver (such as the GBT’s K-band Focal Plane Array) or an interferometer 
(such as the Allen Telescope Array (ATA)). 


Umax (Hz/s) | Sensitivity Loss Factor | Number of Hits in Bin [% of total] | Number of Events in Bin [% of total] 
0.15 1 1838435 [53.1%] 228741 [54.7%] 
0.3 1.4 1216882 [35.1%] 129808 [31.0%] 
0.6 2 376107 [10.9%] 46390 [11.1%] 
1:2 2.8 21999 [0.64%] 7569 [1.81%] 
2.4 5.7 5549 [0.16%] 3232 [0.77%] 
4.8 8 2607 [0.08%] 1514 [0.36%] 
9.6 11.3 816 [0.02%] 441 [0.11%] 
19.2 16 251 [< 0.01%] 134 [0.03%] 
38.4 22.6 139 [< 0.01%] 68 [0.02%] 
76.8 32 172 [< 0.01%] 90 [0.02%] 
153.6 45.3 125 [< 0.01%] 71 [0.02%] 
307.2 64 83 [< 0.01%] 54 [0.01%] 
614.4 90.5 0 [< 0.01%] 0 [< 0.01%] 


Table 3. Maximum drift rates (absolute value), the corresponding multiplier for reduction in sensitivity, and the number of 
hits affected by that reduction. The reason for the sensitivity reduction is described in Section 4. 99% of hits and 97% of events 
that we detected suffer a sensitivity loss of 2 or less. We can also consider the implied population of high-drift hits that we 
did not detect, due to the sensitivity loss. However, the number of high-drift rates is so many orders-of-magnitude below the 
low-drift hits, that we do not expect the number of “missed” hits to meaningfully affect the total number. This implies that 
the sensitivity loss at high-drift rates should not greatly affect the conclusions of our survey. Note that this logic assumes that 
high-drift hits are not preferentially likely to be caused by an ETA. 


Also, as Margot et al. (2021) has mentioned, the turboSETI pipeline does not optimize for signals with high drift rates 
by default, such that signals suffer from drift-smearing across frequency bins, lowering the sensitivity. However, we can 
recover as much sensitivity as is physically possible (in an incoherent data product) by performing stepped frequency 
scrunching. The remaining sensitivity loss factors for this data, and how many signals were affected by those factors, 
are shown in Table 3. While we recommend that other turboSETI users also perform frequency scrunching if their 
searches cover applicably high drift rates, we note that each scrunched data product must be searched separately. This 
means that bright-enough hits will be detected redundantly multiple times, at multiple drift rates, as the hits from each 
scrunched data product are produced independently. This provides a strong motivation to integrate drift-rate-stepped 
scrunching into turboSETI in a more robust way. 

Another way to address high-drift signals is to coherently correct the raw voltage data to an accelerating reference 
frame or fiducial drift rate (e.g., +100 Hz/s) before the creation of the reduced and channelized data products (as 
done in pulsar de-dispersion, e.g., Sobey et al. 2022). For signals at the chosen drift rate, there will be no sensitivity 
loss, and signals that are nearby in drift rate space can be searched for incoherently (with the same sensitivity loss 
considerations as discussed above), broadening the applicability of the technique. However, note that for lower loss 
tolerances at any given drift rate, the problem gets computationally heavier — to minimize loss, it is necessary to 
dedisperse and save the entire array at every drift rate of interest. 

This search, synchronized with planetary transits, is part of a class of search strategies aligned in time with particular 
astronomical events, as discussed in Section 1. For any such search strategy, we expect that the transmitter will likely 
not be transmitting continuously, but instead will only transmit at particular times, potentially with short durations. 
Therefore, searches for these kinds of signals benefit from observing with arrays that can do simultaneous off- and 
on-sources with multiple beams, instead of single-dish single-pixel instruments like the GBT at L-band. 

The short transmission times also make synchronization strategies more energy-efficient for the transmitter, similar 
to the argument made by Gray (2020). Assume that the energy required to transmit for a full exoplanetary orbit is E. 
For the exoplanets in our sample, transmitting only during transit (as seen by a single receiver, in this case Earth), 
costs only 1-10% of E. If the transmitter instead only signalled during the 5 minutes around the mid-transit point 
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(corresponding to the length of one of the observations in this campaign), the energy cost would be 0.1-0.01% of E. 
Another strategy to improve energy efficiency could be to use a transmitter with a very large effective area, leading to 
an extremely narrow beam. In this strategy, the tiny beam size would cause a short flash right at mid-transit as seen 
from Earth. Having such a tiny angular beam size only works well if the observer knows exactly when to expect the 
signal, so it pairs well with transit synchronization. To an order-of-magnitude, using exoplanetary properties consistent 
with our 12 exoplanet sample, we can imagine a continuously-emitting transmitter with an effective area equivalent 
to an exoplanet’s projected area. Before, we assumed the effective area would be the same as that of the GBT, so 
the new effective area is now a billion times larger. This factor propagates, leading to such a small beamsize that 
a transmitted signal would be visible for a few milliseconds at mid-transit, and could be sent with a transmitter six 
orders-of-magnitude less powerful than Arecibo was. This paper does not account for an optimization this extreme, 
but it could be an interesting avenue for more specialized searches in the future. 


7. CONCLUSION 


In this work, we describe the first radio technosignature search that pre-planned observations to synchronize with 
exoplanets during their transits, in a survey of a dozen exoplanets in the Kepler field. Using 6 hours of L-band (1.1- 
1.9GHz) data taken with the BL bankend on the GBT, we performed a SETI search using the narrowband signal 
search code turboSETI. We chose a maximum drift rate of +614.4Hz/s — the first modern radio technosignature 
project to encompass such extreme drift rates — in order to account for the full range of drifts that could be produced 
from known exoplanetary systems. We also chose a low S/N of 5. With these parameters, the algorithms flagged 
~ 2.53 x 10° hits, which were then temporally filtered with an on-off method into 3.4 x 10° events. Many of these 
events could be attributed to GPS satellite downlinks. 

Thirteen citizen scientists volunteered their time to assist the science team with the further filtering of the turboSETI 
events. From this process, the list of signals-of-interest was reduced to a few thousand signals that appeared to be 
either transient or sky-localized. We further removed signals that appeared in multiple targets, or that were identical 
to signals proven to be RFI, reducing the pool further to 479 signals-of-interest. Upon investigating these 479 signals 
in the context of the entire observing session, we determined that only two remained as signals-of-interest: one at 
1749 MHz in Kepler-1332b, and the other from 1040-1438 MHz in Kepler-842b. 

These signals do not rise to the level of even “candidate” technosignature signals because there is no proof that they 
are spatially isolated (and they are consistent with anthropogenic RFI), and so do not warrant follow-up with the 
rigor described in (Sheikh et al. 2021). Reobservation of these targets during transits, with a multibeam instrument 
such as the ATA, will conclude the experiment and, if nothing is found, complete the null result we report here. We 
hope that the “new ground” that we have broken in radio technosignature parameter space will be extended by more 
synchronized SETI searches in the future, across many more instruments and teams. 
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APPENDIX 


A. OBSERVATION LOG 


Table 4. The observation log from this project’s Green Bank Telescope observations. The scan ID number is shown in the 
first column, the target is shown in the second column, the transit synchronization is shown in the third column, and any 
additional notes are shown in the fourth column. Scan IDs are mostly consecutive, but missing ID numbers indicate scans 
which did not take data. Targets prefixed with a “K” indicate a “Kepler” system, with the letter suffix indicating which planet 
in the system was predicted to be transiting during our observing window. The relationship of the observation time to the 
transit time is indicated by letters: I-M indicates an observation which occurred between a transit’s ingress and midpoint, M 
indicates an observation which contained the transit midpoint, and M-E indicates an observation which occurred between a 
transit’s midpoint and egress (blank entries occurred outside of transit). Scan 0051 is missing a node of data (about 1/8th of 
the bandwidth) due to faulty data collection, encompassing frequencies from 1.25-1.30 GHz (entirely within the notch filter). 


Scan ID Target Relation to Transit Notes Scan ID Target Relation to Transit Notes 
0006 30295 non-planet obs. Quasar 0033 K723b LM 
0007 3C295 non-planet obs. Quasar 0034 K1164b M-E 
0009 Boyajian’s Star non-planet obs. Off-source 0035 K723b M-E 
0010 K992b M 0036 K537b M 
0011 K738b I-M 0037 K723b I-M 
0012 K992b M-E 0038 K537b M-E 
0013 K738b I-M 0039 K723b I-M 
0014 K1039b M 0040 K537b M-E 
0015 K1039b M-E 0041 K723b I-M 
0016 K732c M 0042 K723b M 
0017 K1039b M-E 0043 K1332b 
0018 K738b M 0044 K723b M-E 
0019 K738b M-E 0045 K1332b 
0020 K732c M-E 0046 K723b M-E 
0021 K1164b I-M 0048 K723b M-E 
0022 K732c 0049 K446b I-M 
0023 K1164b I-M 0050 K723b M-E 
0024 K732c 0051 K446b I-M missing node 
0025 K1164b I-M 0052 K446b M 
0026 K1053b M 0053 K1222b I-M 
0027 K738b M-E 0054 K1222b M 
0028 K1164b M 0055 K842b I-M 
0029 K1053b M-E 0056 K842b M 
0030 K738b M-E 0057 K723b M-E 
0031 K1332b M 0058 K723b M-E 
0032 K1164b M-E 0059 B0329+54 non-planet obs. Pulsar 


B. MID-RESOLUTION PLOTS FOR THE TWO SIGNALS-OF-INTEREST 
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Figure 10. The signal-of-interest in Kepler-1332b data shown in the high frequency-resolution data (left) and the mid-resolution 
data (right). The regions of data shown are related as indicated by the white-dashed rectangles, with the right plot zoomed-out 
in frequency and zoomed-in in time. The signal-of-interest is actually composed of two short (<1 second) pulses. 
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Figure 11. The signal-of-interest in Kepler-842b data shown in the high frequency-resolution data (left) and the mid-resolution 
data (right). The regions of data shown are related as indicated by the white-dashed rectangles, with the right plot zoomed-out 
in frequency and zoomed-in in time. The signal-of-interest has a distinctive, periodic structure in frequency and is present across 
the bandwidth. This behaviour is characteristic of RFI. 


